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Abstract 

We study the limit of a kinetic evolution equation involving a small 
parameter and perturbed by a smooth random term which also involves 
the small parameter. Generalizing the classical method of perturbed test 
functions, we show the convergence to the solution of a stochastic diffusion 
equation. 
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1 Introduction 

Our aim in this work is to develop new tools to study the limit of kinetic 
equations to fluid models in the presence of randomness. Without noise, this is 
a thoroughly studied field in the literature. Indeed, kinetic models with small 
parameters appear in various situations and it is important to understand the 
limiting equations which are in general much easier to simulate numerically. 
In this article, we consider the following model problem 

+ -a{v) ■ V^r = + -fm^ in x x K, (1) 

with initial condition 

no) - in X V,, (2) 

where i is a linear operator (see ([3|) below) and a random process depending 
on {t,x) G K+ X T'' (see Section [^^ . We will study the behavior in the limit 
£ — > of its solution f^. 

In the deterministic case = 0, such a problem occurs in various physical sit- 
uations: we refer to |DGPOO| and references therein. The unknown f'^{t, x, v) is 
interpreted as a distribution function of particles, having position x and degrees 
of freedom v at time t. The variable v belongs to a measure space (V, /x) where 
H is a. probability measure. The actual velocity is a{v), where a S L°°(F;M''). 
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The operator L expresses the particle interactions. Here, we consider the most 
basic interaction operator, given by 

Lf= I fdfi-f, feL\V,fi). (3) 
Jv 

Note that L is dissipative since 

- f Lf-fdfi=\\Lf\\l.^y,^^, feL^V^f,). (4) 
Jv 

In the absence of randomness, the density = Jy f^d^ converges to the solution 
of the linear parabolic equation (see section [231 for a precise statement): 

dtp- div(A'V/9) = in R+ x T"^, 

where 

K := / a{v) ® a{v)dfi{v) (5) 
Jv 

is assumed to be positive definite. Wc thus have a diffusion limit in the partial 
differential equation (PDE) sense. 

When a random term with the scaling considered here is added to a differen- 
tial equation, it is classical that, at the limit e — > 0, a stochastic differential 
equation with time white noise is obtained. This is also called a diffusion limit 
in the probabilistic language, since the solution of such a stochastic differential 
equation is generally called a diffusion. Such convergence has been proved ini- 
tially by Khasminskii |Has66a[ IHas66b| and then, using the martingale approach 
and perturbed test functions, in the classical article jPSV77| (see also |EK86j . 
|F(;PS07j . |Kus840 . 

The goal of the present article is twofold. First, wc generalize the perturbed 
test function method to the context of a PDE and develop some tools for that. 
We believe that they will be of interest for future articles dealing with more 
complex PDEs. Second, we simultaneously take the diffusion limit in the PDE 
and in the probabilistic sense. This is certainly relevant in a situation where 
a noise with a correlation in time of the same order as a typical length of the 
deterministic mechanism is taken. Our main result states that under some 
assumptions on the random term m, in particular that it satisfies some mixing 
properties, the density p'^ — Jy f'dp converges to the solution of the stochastic 
partial differential equation 

dp = dw{K\7p)dt + po Q^/^dW{t), in R+ x T'', 

where K is as above, VF is a Wiener process in i^(T'') and the covariance 
operator Q can be written in terms of m. As is usual in the context of diffusion 
limit, the stochastic equation involves a Stratonovitch product. 

As already mentioned, we use the concept of solution in the martingale sense. 
This means that the distribution of the process satisfies an equation written in 
terms of the generator (see section [3^ for instance). This generator acts on test 
functions and the perturbed test function method is a clever way to choose the 
test functions such that one can identify the generator of the limiting equation. 
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Instead of expanding the solution of the random PDE /'^ as is done in a Hilbert 
development in the PDE theory, we work on the test functions acting on the 
distributions of the solutions. 

In section [21 we set some notations, describe precisely the random driving term, 
recall the deterministic result and finally state our main result. Section|3]studies 
the kinetic equation for e fixed. In sectional we build the correctors involved in 
the perturbed test function method and identify the limit generator. Finally, in 
section [SJ we prove our result. We first show a uniform bound on the norm 
of the solutions, prove tightness of the distributions of the solutions and pass to 
the limit in the martingale formulation. 

We are not aware of any result on probabilistic diffusion limit using perturbed 
test functions in the context of PDE, but the recent work jdBGllj (in a context 
of nonlinear Schrodinger equations) and |PP03j (where the underlying PDE 
is parabolic and the limit [e — ?> 0] associated to homogenization effects). A 
diffusion limit is obtained for the nonlinear Schrodinger equation in jMar06) . 
jdBDlOj , |DT10| but there the driving noise is one dimensional and the solution 
of the PDE depends continuously on the noise so that in this case an easier 
argument can be used. Eventually, note that a method of perturbed test func- 
tion has also been introduced in the context of viscosity solutions by Evans in 
jEva89| . Actually, in the case m = 0, i.e. for the deterministic version of ([ij, 
the method of |Eva89| allows to obtain the diffusive (in the PDE sense) limit 
[e — > 0] of ([T]) when the velocity set V is finite. 

2 Preliminary and main result 
2.1 Notations 

We work with PDEs on the torus T'', this means that the space variable x S 
[0, 1]'' and periodic boundary conditions arc considered. The variable v belongs 
to a measure space (V, fi) where /i is a probability measure. Wc shall write for 
simplicity ^ instead of L'^(T'^ x V, dx®d^), its scalar product being denoted by 
(•, •). We use the same notation for the scalar product of i^(T''); note that this 
is consistent since ^liy) = 1. Similarly, we denote by the norm (-u, u)^/^, 

whether u G ,„ or L'^iT'^). We use the Sobolev spaces on the torus H^iT''-). 
For 7 G N, they consist of periodic functions which are in L^(T'') as well as their 
derivatives up to order 7. For general 7 > 0, they are easily defined by Fourier 
series for instance. For 7 < 0, H^iT'^) is the dual of H~^{T'^). Classically, 
for 71 > 72, the injection of H^^i^'^) in H'^'^{T''') is compact. We use also 
i°°(T'^) and W^^'="{T'^), the subspace of L'^{T'^) of functions with derivatives 
in L°°(T''). Finally, L"^ {V \ {T'^)) is the space of functions / of w and x such 
that all derivatives with respect to x are in L^(T'^) and the square of the norm 




d 
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is finite. 
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2.2 The driving random term 

The random term has the scaling 

TO^(t, x) = m ^ 

where m is a stationary process on a probabihty space (17, J^, P) and is adapted 
to a filtration (J"t)igR. Note that rrf is adapted to the filtration (J7)tgR, with 
J-,-2t, t € M. 

Our basic assumption is that, considered as a random process with values in a 
space of spatially dependent functions, m is a stationary homogeneous Markov 
process taking values in a subset E of iy^'°°(T'^). We assume that m is stochas- 
tically continuous. Note that m is supposed not to depend on the variable v. 
The law v of m{t) is supposed to be centered: 

Em(t) = / ndv{n) ^ 0. (6) 
J E 

In fact, we also assume that m is uniformly bounded in W^'°"{T'^) so that E is 
included in a ball of W^-'^'°°(T''). We denote by {Pt)t>o a transition semigroup 
on E associated to m and by M its infinitesimal generator. 

As is usual in the context of diffusion limit, we use the notion of solution of the 
martingale problem and need mixing properties on m. We assume that there is 
a subset Dm of Cb{E), the space of bounded continuous functions on E, such 
that, for every ip G Dm, Mip is well defined and 

'ip{m{t)) - / Mip{m{s))ds 
Jo 

is a continuous and integrable martingale. Moreover, we suppose that m is 
ergodic and satisfies some mixing properties in the sense that there exists a 
subspace of Cb{E) such that for any 9 G ^^a/ the Poisson equation 



Mip = e- / e{n)dv{n) (7) 

JE 

has a unique solution ip € Dm satisfying / ipdi/ — 0. When 9 satisfies 



E 



JE 

edv = 0, (8) 

we denote by M^^6 G -Da/ this solution and assume that it is given by: 

/"OO 

M-^e{n) = Pt9{n)dt. 
Jo 

In particular, we suppose that the above integral is well defined. It implies that 



lim Pte{n) =0, Vn G E. (9) 

We need that 3^m contains sufficiently many functions. In particular, we assume 
that for each x G T**, the evaluation function i/jx defined by ipxin) = n{x), n G 
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E, is in ^Af. Also, we assume that, for any f,gGLl the function ipj^g : 
n ^ (/, ng) is in and we define M~^I from E into W^-°°{T'^) by 

(/,Af-i/(7i).9) :=M-V/.<,W, V/,.9ei2_„. (lo) 

We need that M~^I takes values in a ball of W^'°°{T'^) and take C* large enough 
so that 

||n||H'i.-(T<') < a, p/-i/(n)||vvi.=c,(T.) < a, (11) 

for all n £ E. It is natural to require the following compatibility assumption, 
which would follow from continuity properties of M^^: 

M-^^W = M^^I{n){x), yneE,x e T^. (12) 

Note that by -tf^j^g and '0x satisfy the centering condition ([5]). Note also 
that, by ([TU)) and p^ . we have, taking <? = 1, 

/" f{x)Ar^iP^{n)dx = M~^iPf^x{n). (13) 

Eventually, we will also assume that for any f,g £ ,„ and for x G T"*, the 
functions 

^f,g : n^if,nAr'l{n)g), M-^^f,i, M^Vx 

are in ^m- 

To describe the limit equation, we remark that since m(0) has law v, 

^y{n)M'^^P,{n)dv{n) = -E (^j,(m(0))A'/-V=.(m(0))) (14) 

ijy{m{0))J^ Ptij,{m{0))d?j (15) 

E i^iPy{m{0)) ij,{m{t))d?j (16) 

E (m{0){y) [ 77i{t){x)dt) , (17) 



where we have used the Markov property in the identity ([T5]) - (fTO|) . We define 
k e L°°(T'^ X T'') by the formula 



k{x, y) = E / m{0){y)m{t){x)dt, x, y e T''. 
Let F e L°°(T'^) be the trace 

F{x) = k{x, x)^E m{0){x)m{t){x)dt, x e T'^. 
Note that, m being stationary, 

k{x, y) = E m{0){y)m(t){x)dt^ ^ ^ (/ m{'^){v)'m'{t){x)dt 

= E ■m{0){y)m{t)(x)d?j + ^(^J m{-t){y)m{Q){x)dt 

= E^m(0)(y)y m{t){x)d?j +E(^m{0){x) m{t){y)dt^ , (18) 
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so that k is symmetric. Let Q be the Hnear operator on L^(T'^) associated to 
the kernel k: 



Qf{x)^ / k{x,y)f{y)dy. 

Lemma 1. The operator Q is self-adjoint, compact and non-negative: {Qf, f) > 
for all f £ L^{T^). 

Proof: Q is self-adjoint and compact since k is symmetric and bounded. To 
prove that {Qf, f) > 0, we will need the following fact: if £ satisfies ([5]), 
then 

(19) 



lim E\PTip(m(0))\ = 0. 

T-i>+oo 



E|Pt^(™(0))| = / \PTi^{n)\diy{n) and |Pt^H| < MciE), 



Indeed 



whence (|19p by the mixing property ^ and by the dominated convergence 
Theorem. In particular, if ijj G satisfies (|5]) and if, furthermore, M~^ijj e 
3^M, then 

dt 



hence 



PMmiQ)) 



|PtM"^V'(to(0))| ^ 



(20) 



when T — > +oo. For simplicity, let us denote hy ipf the function V'/.i- By (|13p . 
([T4|. ([T8l) and (EOl), we have 



(Q/,/) = -2E [V'/(™(0))Af-V/(m(0))] 



2E 



V;/(m(0)) / Pt7/;/(TO(0))dt 







0(1), 



when T — > +cx3. On the other hand, for T > 0, we compute 

rr, 2 



-E 



T 







'ip}{m{t))dt 



^0 
T /.t 



2 rT rt 
T 

2 
T 
2 

2^ Jo ^0 
2 
f 



(21) 



E[i/'/(m(t)))i/'/(m(T))]dtdT (22) 

E[?/>/(m(t))?/'/(m(T))]dtdT (23) 
E[V'/(m(t-r))V'/(TO(0))]dtdT (24) 
E[V'/(TO(r))?/'/(m(0))]dtdr 
(T - T)E[V'/(m(T))?A/(m(0))]dT 



"'0 



= 2 E[2jjf{m{T))^jf{m{0))] + rr 



V'/(™(0)) / Pt^/(m(0)) 



2E 



(25) 
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where we have use the homogeneity of m{t) in ([23])- ([Ml)- The remainder 
satisfies 



1 '•^ 



rr - -2E (^^/>y(m(0))- J rPri^ tim{Q))dT 
Since tpf g ^m, -PrV'/ = j^PrM^^^if. this gives 

TT = -2E I V/(to(0)) 
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By (Uni), we obtain rr = o(l). By dH]), (Q/, /) is the hmit of the left-hand side 
of ((22|) . which is non- negative, hence {Qf,f) > 0. ■ 

As a resuh of Lemma[TJ we can define the square root Q^^^. Note that Q^^^ is 
Hilbert-Schmidt on L'^{T'^) and that, denoting by WQ^^^Wj^^ its Hilbcrt-Schmidt 
norm, we have 

||Qi/2|||,^ ^ Tr Q = / k{x,x)dx. 

We will not analyze here in detail which kind of processes satisfies our assump- 
tions. The requirement (jlip that m and M~^m are a.s. bounded in W°°{T'^) 
are quite strong. An example of process we may consider is 



m{t) 



with rij e W^'°°{T'^), 



j6N 



where the processes (mj)jgN are independent real valued centered stationary, 
satisfying the bound 

\m^{t)\ < C, a.s., t e R, 

for a given C > 0. We are then reduced to analysis on a product space. The 
invariant measure of m is then easily constructed from the invariant measures 
of the mj 's. Also, the Poisson equation can be solved provided each Poisson 
equation associated to rrij can be solved. This can easily be seen by working 
first on functions ip depending only on a finite number of j. 
The precise description of the sets Dm and ^^v/ depends on the specific pro- 
cesses TTij, j G N. For instance, if mj are Poisson processes taking values in 
finite sets Sj, then Di\j and can be taken as the set of bounded functions on 
riiGW^j- More general Poisson processes could be considered (see |FGPS07| ). 

Actually, the hypothesis (fTT|) can be slightly relaxed. The boundedness assump- 
tion is used two times. First, in the proof of (pO)) and (|3T|) . but there it would be 
sufficient to know that m has finite exponential moments. It is used in a more 
essential way in Proposition 1101 There, we need that the square of the norm 
of TO and M^^m have some exponential moments. However, (under suitable 
assumptions on the variance of the processes for example), we may consider 
driving random terms given by Gaussian processes, or more generally diffusion 
processes. 
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2.3 The deterministic equation 

There are also some structure hypotheses on the first and second moments of 
fi: we assume 

/ a{v)dn{v) = 0, (26) 
Jv 

and suppose that the following symmetric matrix is definite positive: 

K -.^ a{v) (g> a{v)dfi{v) > 0. (27) 
Jv 

An example of {V,fi,a) satisfying the hypotheses above is given by V = 
(the unit sphere of M.'^) with ji ~ d — 1-dimensional Hausdorff measure and 
a{v) = V. 

In the deterministic case m = 0, the limit problem when e — > is a diffusion 
equation, as asserted in the following theorem. 

Theorem 2 (Diffusion Limit in the deterministic case). Suppose m = 0. As- 
sume that (/q) is bounded in „ and that 

Jv 

Assume i26]) - [Wl\ ). Then the density p'^ := Jy f^dp converges in weak-L^ ^ to 
the solution p to the dijfusion equation 

dtp - div(A'Vp) = M+ X T"*, 

with initial condition: p(0) = po in T^. 

This result is a contained is [DGPOO] where a more general diffusive limit is 
analyzed. Note that, actually, strong convergence of (p^) can be proved by 
using compensated compactness, see |DGP00j also. 

2.4 Main result 

In our context, the limit of the Problem ([I])-© is a stochastic diffusion equation. 

Theorem 3 (Diffusion Limit in the stochastic case). Assume that (/q) is 
bounded in and that 

Jv 

Assume Then, for all r] > 0, the density p" := Jy fdp 

converges in law on C{[0,T]; H^'^) to the solution p to the stochastic diffusion 
equation: 

dp = div{K\/p)dt + ^Fp + pQ^^^dW{t), in M+ x T'', (28) 

with initial condition: p{0) = pQ in T'^. In (j28p . W is a cylindrical Wiener 
process on L^(T''). 
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It is not difficult to see that formally, (pS)) is the Ito form of the Stratonovitch 
equation 

dp = Aw{KVp)dt + po Q^''^dW[t), in R+ x T''. (29) 

Theorem 13] remains true in the slightly more general situation where the coeffi- 
cient in factor of the noise in ([T]) is in the form icr(/)m^ with 

ct(/) = o{p) + /, p := / fdp, 
Jv 

where a is a. smooth, sublinear function. 

3 Resolution of the kinetic Cauchy Problem 
3.1 Pathwise solutions 

Problem ([I])-© is linear and solved for instance as follows. Let A a{v) ■ 
denote the unbounded, skew-adjoint operator on ^ with domain 

D{A) :={/ eLl,;a{v)-V.J eLl,}. 

Since A is closed and densely defined, by the Hille-Yosida Theorem |CH98j . it 
defines a unitary group e*"^ on ^. 

Theorem 4. Assume Then, for any /q G ^ and T > 0, there exists 

a unique solution F-a.s. in C([0,T];L^^) of ([Ip-(0) on [0,T], in the sense 
that, 

P-a.s., for all t G [0,r]. Besides, if f^ 6 L'^{V; H^{T'^)), then, P-a.s. f 6 
C^i[0,T];Ll^)nC{[0,T];L^{V;H\T'^))). 

The proof of this result is not difficult and left to the reader. The last statement 
is easily obtained since A commutes with derivatives with respect to x. 

Energy estimates can be obtained. Indeed, for smooth integrable solutions 
to (dl-©, we have the a priori estimate 

jjrmh-^iLr,n = -J(a(«)-vr,r) + J(/w,r) 

£ 

By and ([TT|). this gives the bound 

\\rmh + - \\Lf%s)\\i.ds<m\\i. + ^ wrmh ds, 

£ Jo £ Jo 

hence, by Gronwall's Lemma, the following bound (depending on s): 

Iir(t)|li. <e^*||/o1li.. (30) 

Similarly, we have 

\\nt)\\hHr)<e'^'m\\hiHn- (31) 
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It is sufficient to assume G {V ; {T'')) {resp. G (V ; H'^ (T^))) to 
prove ([50)1 {resp. (PT|) ). By density, the inequality holds true for /q G L^ ,, 
(resp. G L^(F;iJ^(T''))). In particular, ||/''(t)||L2 is uniformly bounded in 
w G r! if G „ and 11^(^)11^2(^1) also if G L'^iV; H^T^)). 

3.2 Generator 

The process is not Markov but the couple {f^,m^) is. Its infinitesimal gen- 
erator is given by : 

^'ifi = -^A^'P + \^L*V, (32) 

with 

^aMI, n) = -{Af, Dipif, n)) + {fn, Dip{f, n)), 

^lMI. n) = {Lf, D^if, n)) + A/^(/, n). 

These are differential operators with respect to the variables f € L"^ n ^ E. 
Here and in the following, D denotes differentiation with respect to / and we 
identify the differential with the gradient. For a function on we also 

use the second differential D^tp of a function ip, it is a bilinear form and we 
sometimes identify it with a bilinear operator on by the formula: 

D\{f)-{h,k)^{D\{f)h,k). 

Let us define a set of test functions for the martingale problem associated to 
the generator 

Definition 5. We say that 'J is a good test function if 

• : L'^{V;H\T'^)) x E ^ R, {f,m) ^ 'i>{f,m) is differentiable with 
respect to f 

• (/, to) ^ D'^{f,m) is continuous from {V \ {T'^)) x E to „ and 
maps bounded sets onto bounded sets 

• {f,m) 1^ M'i'{f,m) is continuous from L'^ {V ; {T'^)) x E to M and maps 
bounded sets onto bounded sets of M 

. for any f eL^{V;H\T^)), ■^{f,-)^DM. 

We have the following result. 

Proposition 6. Let 'if be a good test function. Let /q G L'^{V; H^{T'^)) and let 
f^ be the solution to Problem ([2P-(0). Then 

MUt) ■.^^{rit),m'it))- f J^^-^{r{s),m'{s))ds 

Jo 

is a continuous and integrable (J-f) martingale with quadratic variation 

(M|,M|)(t)= / (if^|*p-2*^^*)(/-^(s),TO^(s))ds. (33) 
Jo 
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Proof: Let s,t > and let s = < • • • < t„ = t be a subdivision of [s, t] such 
that maxi — ti \ = S. We have for any J-"J measurable and bounded g 

E(^(^*(.r(i),™ni)) - vI/(r(s),TO^'(s)))5' 

With 



^ = ^e( ( - *(r(<,),m^(t,+i)) 



and 



n-l 



We write 



with 



A = E 



as{s)ds ]g], 



1=1 

Since /q G ^(T'')), we deduce from (pij) and the assumption on ^I* that 

ai is uniformly integrable with respect to (s,a;). Also is almost surely contin- 
uous and TO^ is stochastically continuous. It follows that D'i {f^ (s) , m'^ (ti+i)) — 
{f^ (s) , m'^ {s)) converges to in probability when 6 goes to zero for any s. 
By uniform integrability, we deduce that A converges to 0. Similarly, we have 

t,+i 



{U),m'{a))-M^{r{a),m'ia))da g 



and, by the same argument, B converges to zero when S goes to zero. The result 
follows : M|, is a continuous martingale. Since 5' is a good test function and /q £ 
L^{V; H^iT'^)), it follows from ^ and the bound (HH that t ^ ^'(/^(i), "i^W) 
and t 1-^ ^^'^{f^{t),m'^{t)) are a.s. bounded. The expression for the 
quadratic variation can then either be computed by expanding 



E|M|,(t)|2 =E 



n— 1 



^'^{f'{s),m^{s))ds 



where = < ■ • • < i„ = t is an arbitrary subdivision of [0, t] with step (5 J, 0, 
or, quite similarly, by proceeding as in Appendix 6.9.1 in [FGPS07| . ■ 
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4 The limit generator 



To prove the convergence of (p^), we use the method of the perturbed test- 
function }PSV77| . The method of |PSV77j has two steps: first construct a 
corrector <~p^ to (/? so that ££'^'~p^ is controlled, then, in a second step, use this 
with particular test- functions to show the tightness of (p^). In the first step, we 
are led to identify the limit generator acting on (p. 



4.1 Correctors 

In this section, we try to understand the limit equation at e — !■ 0. To that pur- 
pose, we investigate the limit of the generator by the method of perturbed 
test-function. 

We restrict our study to smooth test functions and introduce the following class 
of functions. Let ip G C^{L'^ y). We say that (p is regularizing and subquadratic 
if there exists a constant C^p > such that 



f l^(/)l<c^(i + ll/IU=) 



2 



||A™z?<^(/)IU. <c^(i + ||/IU.), 

■ iA"^'h,A"^^k)\ < CJh\\L2\\k\\L2, 



(34) 

|DV(/) • iA"^'h,A"^-k,A"'H)\ < CJh\\L2\\k\\L4l\\L-, 



for all f,h,k,l e L^^„, for aU m,mi G {O,-- - ,3}, i = {1,2,3}. Note that 
regularizing and subquadratic functions define good test functions (depending 
on / only). 

Given ip regularizing and subquadratic, we want to construct p>i, p2 good test 
functions, such that 

ifV"(/, ") = ^<^(/, n) + 0{e), ^'^^ + epi+ £^2- 

The limit generator ^ is to be determined. By the decomposition ((32)) . this is 
equivalent to the system of equations 

^L*(p = 0, (35a) 

^A*<^ + ^L*<^i -0, (35b) 

^A*^i + - i?<^(/, n), (35c) 

=S^A,¥'2 = 0(\). (35d) 



4.1.1 Order e ^ 

Equation p5a[) constrains pi to depends on p = f ~ Jy fdfi uniquely: 

p{f) = p{p), p:= f fdfi, (36) 
Jv 

and imposes that the limit generator ^ acts on p{f) uniquely, as expected in 
the diffusive limit, in which we obtain an equation on the unknown Jy fdji. 
Indeed, since p} is independent on n, p5ap reads 

{Lf,Dp{f))^Q. (37) 
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Let {g{t, f))t>o denote the flow of L on L'^{V,fi): 

j^g{tJ)^LgitJ), 3(0,/) = /. (38) 
An explicit expression for g is 

gitj)^p + e-'if~p), p^ f f{v)dp{v). 

Jv 

In particular, g(t, f) ^ p exponentially fast in L'^{V, i-i) when t +00. By 
equation (j37p is equivalent to 

^U) = ^{g{tj)), vtGR, 

i.e. (j36p by letting t — > +00. 



4.1.2 Order 

Let us now solve the second equation (j35bp . To that purpose, we need to invert 
^L*- Let us work formally in a first step to derive a solution. Assume that 
m{t,n) is a Markov process with generator M, let g be defined by (p8|) and 
consider the Markov process {g{t, f),m{t,n)). Its generator is precisely ^l*- 
Denote by {Qt)t>o its transition semigroup. Since both g and m satisfy mixing 
properties, the couple (g, m) also. In particular, we have 



QMf,n) ^ (V')(/) := / ^{f,n)d,y{n), (39) 



E 



and it is expected that, under the necessary condition {^Atf) = 0, a solution 
to (|35bp is given by 



/>oo 

Jo 



Let us now compute ^a*'P- By we have for h G {h,Dip{f)) — 

{h., Dip{p)), where as above the upper bar denotes the average with respect to 
V and p :~ f. Hence 

^AMf,n) = -'(Af,Dip{p)) + ipn,Dip{p)). 

Since the first moments of a{v) and m{t) vanish, we have 



Ap = and / {pn, DLp{p))dv{n) = 0, 
J E 

and the cancellation condition {J^a*^) = is satisfied. We then write 

/•oo 

<pi(/,n)= / Qt^A*f{f,n)dt 

J 

E{^AM9itJ),Mt,n)))dt. 





00 
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Note that g is deterministic and g = p, so that 



(pi(/,n)= / ~iAg{tJ),Difip))+Eiipmit,n),D^ip)))dt 



(Agit, f),Dip{p))dt - (pAr'lin), D^{p)). 



Furthermore, regarding the term Ag{t, /), we have 

j^Ag{t, f) = A^^git, /) = ALg{t, f) = Agit, f) - Agit, /). 



Since Af = 0, we obtain —Agit, f) = —Agit, /), i.e. 

at 



Agit,f) = e-'Af. 

It follows that 

Mf,n) = ~iAf,D^ip)) - ipM-Hin),D^ip)). 
By this is also equivalent to 

Mf,n) = -iAf,D^if)) - ifM-'lin),D^if)). (40) 
This computation is formal but it is now easy to define ipi by (j40[l and to check 



that it satisfies (j35bp . It is also clear that is a good test function. 

Proposition 7 (First corrector). Let ip G C^iL'^^,) be regularizing and sub- 
quadratic according to ^3^. Assume that (p satisfy Ii36\} . Then i35b]) has a 
solution pi G C^iL^ ^ x i?) given by 

p,if,n) = -iAf,Dpif)) - ifM-'lin),Dpif)), (41) 
for all f £ L/^ J,, n G E. Moreover pi is a good test function. 

4.1.3 Order e° 

Let us now analyze Equation pScp . Setting p — f, it gives 



^pip) = {^A*Pi)ip) = / .^A*Mp,n)diyin). (42) 

JE 

We have 

^aMI, n) = i-Af + fn, D^if, n)) 

and 

Mf,n) = i-Af - fM-^Iin),D^if)) = -iAf,Dpip)) - ipAr'lin), Dpip)) 

p{if,n) +plif, n). 

By (|42p. the limit generator is therefore the sum of two terms: 

^pip)^j^Mp) + -^Mp)- 
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The first term ^^ip[p) corresponds to the deterministic part of the equation. 
We compute, for h £ 

(/i, = -(Ah,D^{p))-D\{p) ■ (Af, h). 

In particular, evaluating at / = p we have 

{h,D4{p,n)) ^ ~(Ah,D^{p)) 

since Ap = 0. Taking then h — —Ap + pn and using once again the cancellation 
property Ap = 0, we obtain 

^Mp)^ I (A^P,D^{p))d,y{n), 



I.e. 



J^i^{p)^{A^p,D^{p)). (43) 

The second part .Sf* corresponds to the random part of the equation: since 
Tp = Q, 

^Mp)= I [pn,D^l{p,n))dy[n), ^\{f ,n) ^ ~{pM-^ I{n), D^{p)). (44) 

JE 

Now that = {^A*Vi) has been identified, we go on with the resolution of 
3c|) . At least formally at a first stage, we can set 



To the decomposition (pi — {p\ + ip*i then corresponds a similar decomposition 

for ip2. Since ^\{n) is linear with respect to n, the term 

^A*^l{f,n) := {-Af + fn,D^l{f,n)) 

can be decomposed into two parts: one that is linear with respect to n, the 
second that is quadratic in n. The first (linear) part does not contribute to 
since m{t) is centered: Em(t) = 0. Let us thus compute the remaining part 

(l{f,n) := {fn,D(pl{f,n)). 

Since Lp*i{f,n) depends on p only, we have q{f,n) — {pn, Dip*i{p,n)). Since, 
along the flow of i, the density g{t, f) ~ p is constant, we obtain 

^*{f,n) ^- 1^ Pt i^j^q{p,n)dy{n) - g(p,.)| {n)dt. 

In particular, from the expression (j44[) for if \ and the fact that (p is subquadratic 
and regularizing, it follows that G C{L'^ ^ x i5) is a good test function and 
satisfies 

|^^(/,r^)| <C(l + ||/||i.), (45) 
\.^A.V2U.n)\<C{l + \\f\\l.), (46) 
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for all f € ^, n € E, where C is a constant depending on the constant C* 

in and on the constant Cip. Similarly, j2fA*(p\{f, n) is the sum of one term 
independent on n and one term linear with respect to n. This latter does not 
contribute to Lp\ by the centering condition 'Em{t) = 0. We explicitly compute 
the first term: 

{~Af,D^\{f,n)) = {W],D^{p))+D\{p) ■ (Af,Af). 

We have already proved (cf Section 14.1.2^ that, along the flow g{t,f) of L, we 
have Ag{t, /) = e~*Af. Similarly, we have 

A^git, f)-^= e-*(IV - 
By integrating the exponential e~* with respect to t, it follows that 

4if,n) = (^^^p,Dipip)) + ^D'ip{p)(Af,Af). 

In particular, G „ x iJ) is a good test function and satisfies (pSjl -ipS ]) . 

By introducing Ap = div(/\ Vp), where K is given by ([S]), to identify Jfj in (^5]). 
and by developing the expression (|44|) for we obtain the following result. 

Proposition 8 (Second corrector). Let (p £ C'^(L^ „) be regularizing and sub- 
quadratic according to I^J). Assume Ii36]) and ill}) . h26\) . \2T^ . Let A 
denote the unbounded operator defined by 

Ap = div(A'Vp), D{A) = H'^{T^) c L^{T'^). 

Then mdjl is satisfied for ^ defined by: Mi' e C^iL^iT'^)), 

^^{p) = {Ap,Di^ip)) 

-[ {{pnM-^I{n),Dip{p))+D^ip{p)-{pM-^L{n),pn)}di^{n), (47) 
Je 

and a corrector ip2 € C (i^ ^, x E) which is a good test function and satisfies 

\Mf,n)\<C{l + \\f\\l.), 
\J^A*^2{f,n)\<C{l + \\f\\l.), 

for all / e n G E , where C is a constant depending on the constant C* in 

^11\) and on the constant . 

4.2 Limit equation 

We will show here that ^» is the generator of the semi-group associated to a 
diffusion process on L^(T'^). Then (PH) is a form of corresponding to the 
Stratonovitch formulation of the corresponding stochastic differential equation. 
Actually, we use the expanded form of (j47|) (which corresponds to a stochastic 
differential equation in Ito form) to identify more precisely the limit generator 
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S£ . The notations for F, fc, Q are those introduced in Section 12.21 We have 
first: 

/•OO 

(^pnM-^I{n),Dif{p))dv(n) \ {pm{0)m{t), Dip{p))dt 

E Jo 

= \e f {pm{Q)m{t),D^{p))dt 
^^{pF,Dipip)), 

where 

F{x) :=E / m{0){x)m{t){x)dt = k{x,x). 
Jr 

To recognize the part containing D^ip, we identify D^cp with its Hessian and 
first assume that it is associated to a kernel 4>. Then, we write: 

D'^ip{p) ■ {pM~^n, pn)dv{n) 

E 

OO 



= E / D^ip{p) ■ {pm{t), pm{0))dt 
Jo 

= ^Ej^D'ipip) ■ {pm{t),pm{0))dt 

= ^eJ {D\{p){pm{t)),pm{Q))dt 

_ 1 

^ 2 

_ 1 

~ 2 



$(a-, y) p{x)m{t){x) p{y)m{Q){y)dxdydt 



2/)fc(a;, y)p{x)p{y)dxdy. 

fd J Yd 



Denote by q the kernel of Q^/^, then 



Kx^y) = / qix,z)q{y,z)dz, 

Jfd 



which gives 



I?V(p) • {pM~^n,pn)dv{n) 

E 

1 



2 J'J'd 



$(a;, j/)g(a;, z)p{x)p{y)dxdydz 

Yd Jjd 

= \Tmce[{pQ'^')DMp){pQ'^'n (48) 



By approximation, this formula holds for all function ip. We conclude that 
^ is the generator associated to the stochastic PDE 

dp = div{KVp)dt + ^Fpdt + pQ^/^dW{t), (49) 

where is a cylindrical Wiener process. 
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4.3 Summary 

By Proposition [BJ Proposition [7] and Proposition [51 wc deduce: 

Corollary 9. Let ip G C^{L'^^) be a regularizing and subquadratic function 
satisfying i36\) . There exist two good test functions ipi, ip2 such that, defining 



ip" = if + eifi + e'^ip2, 

|^i(/,n)|<C(l + ||/||i.), 

|^2(/,n)| <C(l+||/||i.), (50) 

|^-^-(/, n) - if^(/, n)\<Ce{l + ||/||i.) , 

for all f ^ L^. ^, n G E, where C is a constant depending on the constant C* in 
ill]) and Cip. Moreover 

M'{t) ip%r{t),m'{t)) - f ^'cp'{ris),m'is))ds, t > 0, 

is a continuous integrable martingale for the filtration {J-f) generated by with 
quadratic variation 

{M',M')it)^ f {iM\ipi\^ -2ipiMip,){r{s),m'{s))+r,it)}ds, (51) 
Jo 

where 

\re{t)\<Ce [\l + \\f^t)\\l.)ds, (52) 
Jo 

for a constant C depending on C* and C^p. Finally, for < si < • • ■ < s„ < 
s<t andi^ eCbULl,^)"), 



E l^{p-(t))~ip{p%s))- / ^ip{p%a))da]^P{p%si),...,p%s„)) 



<Ce{l+ sup E\\f%t)\\l. , (53) 

with another constant C depending on the constant C, in ill]) , on Cip and on 
the supremum ofip. 



Proof: Everything has already been proved except for (ISTI) and the last state- 
ment (l53l). For this latter, it suffices to write: 



ip{p%t))-ip{p-{s))- f' ^ip{p^{a))da 



Ar{t) - M^{s) - sMp'it)) - e^Mp'it)) + epi{p%s)) + e^^ip^s)) 
f f^^(p^(a))-if>^(p^(a))V<^- 
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Then, wc multiply by 'ip{p'^(si), . . . , p'^(sn)), take the expectation and use the 
bounds (jSO]) to conclude. Furthermore, 

M\ipf - 2ipMip = (54) 

if (y9 A^(y9 is a linear first order operator in ip. Applying (j54p to 

Mifiif, n) = -^aMJ, n) + ^{L,f, Dip{f, n)) 

gives 

^e|(^E|2 _ 2y3=ifV^ = M|v?ip - 2v5iM(pi +re. 
By ([5n|) . the remainder satisfies ([5^ . ■ 

5 Diffusive limit 

Our aim now is to prove the convergence in law of p'^ = f^d^L to p, solution to 
(P^ . or equivalently Equation To that purpose, we use again the perturbed 
test function method to get a bound on the solutions in L^^ then we prove that 
p^ is tight in C{[Q,T];H~'^), rj > 0. We follow (and adapt to our context) 
the method in |FGPS07| . paragraph 6.3.5. In particular, we use Kolmogorov 
criterion to get tightness in section [221 an alternative method would be to use 
Aldous' criterion for tightness (e.g. Theorem 4.5 in |JS03j ). 

5.1 Bound in L^^ 

Proposition 10 (Uniform ,„ bound). Assume Ulp . Then, for all T > Q, 
p > 1, we have 

E sup wrmi, < c 

te[o.T] 

where the constant C > depends on T, on p, on ||a||^oc(y-), on the constant 
C* in i fTT]) and supg^Q ll/o 11^2 only. 

Proof: Fix p > 2. Let us write a{e,t) < b{e,t) if there exists a constant C 
depending on T, on p, on ||a||Loo(y) and on the constant C* in (fTTj) only such 
that a{e,t) < Cb{£,t) for all t e [0,T]. Set ip{f) := 111/11^2. We want to apply 
Corollary [3] to ip. This requires some care since ip is actually a function of / and 
not of p. Thus, we first seek for one corrector ipi £ C'^{L'^ v ^ ^) such that, for 
the modified test-function 

ip^ ■.= ip + Eipi, 

the term .if^(p^(/^ m^) can be accurately controlled: for / e L^{V; H^iT"^)), 
n ^ E, we compute 

J^>^(/, n) = e-^^lMI) + £~\^A.^ + ^L.^i){f, n) + ^A*Mf, n). (55) 

Since Mip{f,n) — {(p being independent on n), and since Dip{f,n) — /, the 
first term in ((55)) is 

e-2^i,^(/) = -ll!L/||i., (56) 
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which has a favorable sign. Since A is skew-symmetric, J^'A*(p{f,'n) ~ {fn,f). 
This term is difficult to control and we choose (pi to compensate it. We set 

so that Mipi = —{fn, f) and the second term in ([55]) is 

e-\^A*V + ^L*Vi) - -{Lf,D^AJ,n)) - --{LfJM-^I{n)). 

£ e 

By ((TT|) . we obtain 

The remainder JfA*(pi satisfies the following bounds 

^A*Mf, n) = -{Af, fM-H{n)) + (/n, fM-'l{n)) 
= ^{f,AAr'l{n)) + {fn, fM-'l{n)) 

By (ini), (ESI), dSIl), we obtain 
Set 

M^{t) :=^^(r(^),™^(^))-¥>^(/o^m-(0))- ri^V(r(s),m^(s))ds. 

Jo 

By definition of (ys, iy9^ and R-P , we have 

= l\\mli^~eiMrit),m%t)))-Mf^,m^O))) 

Jo 

By ([55)1 and the estimate 

i^i(r(t),m^^W)i<iirwiii.^„, (59) 

we deduce the bound 

wrmh <m\\h +e\\rmi. +f\\rmh ds + snp \M%t)\. 

Jo ' *e[o,T] 

For e small enough, it follows that 

wrmh <m\\h + fwrmh ds+ sup ia^wi, 

20 



and, by Gronwall Lemma, 

wrmh <m\\h + sup \M^t)\. m 

te[o,T] 

On the other hand, similarly to (|5T]) . we have 

{M',M'){t)^ f {M\ipi\^ ~2ipiMipi)iris),m%s))ds. 
Jo 

(Note that there is no remaining terms here since = 0, cf. the proof of ((5T|) 
in Corollary [S]) In particular, by (|59l) and the similar estimate for M(pi, we 
have 

\{M^Amt)\< f\\ns)\\i, ds. 

Jo 

Since M' is a martingale with EM^(t) = 0, Burkholdcr-Davis-Gundy inequality 
gives 

E[ sup \M^t)n<CpE\{M-,Rn{T)\P/' < r E\\r{s)\\f, ds. (61) 
te[o,T] Jo 

By nnm%^^ < nmli^+nsv^Ptmr] iM^im. Hence, by m, 
nr{T)\\% <E\\f^\\i^, +fnr{s)\\%ds. 

Jo 

By Gronwall Lemma, we obtain E||/^(T)||f2 < EH/o^Hf^ . This actually 
holds true for any t G [0,T]. Thus, using (|6T|) and then ((60|) gives finally 

E[sup,,[o,T]iirwiii^ ]<nm% •■ 



5.2 Tightness 

Proposition 11 (Tightness). Let T > 0, ij > 0. Assume (0), ^J\), 
and assume that (/q ) is bounded in . Then {p^) is tight in C([0, T]; iJ~''(T )) 



Proof: Let 'p(p) = Pi or, more precisely (since the perturbed test-function 
method has been developed for real- valued, regularizing functions), define the 
test-function ipj as follows. Let {pj',j > 1} be a complete orthonormal system 
in L^Cfrf)^ let 7 > max{3,4 and let 

J = (Id-A,)~^/^ 

where Id is the identity on L^(T''). Note that ||/5||h-t(t<i) = II^'^^pIIl^ and that 
is Hilbert-Schmidt on L'^{T'^) since -f > d. We set (pj{p) = {J''p,pj). It is clear 
that tpj is subquadratic (it is linear) and regularizing as in p4p (the operator 
V^J'^ is of order < 0). Let be the correction of ipj given by Corollary |9l Let 

M^it)=cp-^{r{t),m%t))~cp-^{f^,m^{0)) - f ^'v^{r(s),m-{s))ds 

Jo 
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and 



Then 

By the estimates ([50]) on the correctors of (pj and the L^-bounds of Proposi- 
tion [TUl we deduce that 

E[ sup \vApHt))-0',m<e, (62) 

*6[0,T] 

where we write a{j, e) < b{j, e) if there exists a constant C depending on 
on T, on the constant C» in (jlip and on sup^^Q ||/o Hl^, but not on s 
and j such that a(j, e) < Cb{j,e). Note also that, by ([50)) . EsupjgjQ 16*^(^)1 ^ 
1, hence 

is a.s. well defined for all t G [0,r] in H '(T'*) since the sum is convergent in 
L2(T''). By dill), we obtain 

E[ sup \\p'(t)~e-mH-,(T^)]<e. (63) 
te[o,T] 

Let ?7 > 0. Let 

w{p,5):^ sup ||p(t) - p(s)||^-„(T<i) 

\t-s\<S 

denote the modulus of continuity of a function p e C([0, T]; i7^''(T'^)). Since 
the injection L'^{T^) C iJ"''(T'^) is compact, the set 

Kn = \peC{[0,T]-H-\T''))- sup ||p(i)||L^ < i?; z«(p, 5) < £(5) I , 
{ *e[o,T] J 

where i? > and e{5) when (5 0, is compact in C([0, T]; H-''\T'^)) (Ascoh's 
Theorem). By Prokhorov's Theorem, the tightness of {p') will follow if we prove 
that, for all a > 0, there exists i? > 0, such that 

P( sup \\p'{t)\\L^ >R)<a, (64) 
te[o,T] 

and 

lim limsupP(z/;(p^, (5) > a) = 0. (65) 
The estimate (|64| follows from the L^-bound of Proposition [10] by the estimate 

P( sup \\p%t)h2 >R)< sup \\p%t)\\L2]. 
t£[0,T] -ft te[o,T] 

Similarly, we will deduce ((65)) from a bound on Eiy(p^,(5) for S > 0. Actually, 
by the L^-bound of Proposition [TU] and by interpolation, we have 

E sup ||p(f) -p(s)|| ,b, <E sup \\pit) - pis)\\'' i 
\t-s\<s ^ ' \t-s\<s 
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for a certain cr > if 77' > jy^ > 0. Therefore it is indeed sufficient to work with 
r] ~ "f. Besides, by (|63p . it is sufficient to obtain an estimate on the increments 
of 9'^ . By definition 



J s 

for < s < t < T. By the L^-bound of Proposition [TUl and by ([50]) . we have 



E 



<\t~s\\ 



By Burkholder-Davis-Gundy inequahty, 

where {Mj,Mj) is the quadratic variation of AfJ. By ([ST|) and the L^-bound of 
Proposition [TUl we obtain 

E|MJ(t)-Af;(s)|4< 
Finafiy, we have E|6'|(t) - 6'|(s)|'' < \t - sp, and tlius 

It foUows (by the Kolmogorov's criterion) that, for a < 1/2, 



By the embedding 



W^"^''(0,r;i7-'^(T'^)) C C"-^([0,T];ff-''(T'^)), < a - i, 

we obtain Ew{d^ , 6) < (5^ for a certain positive fi. This concludes the proof of 
the proposition. ■ 



5.3 Convergence 

We conclude here the proof of Theorem [31 Fix 77 > 0. By Proposition [TTl 
there is a subsequence still denoted by (p^) and a probability measure P on 
C{[0,T];H-''i{T^)) such that 

weakly on C([0,T]; ff-''(T'')) 

where P'^ is the law of p'^ . We then show that P is a solution of the martingale 
problem, with a set of test functions specified below, associated to the limit 
equation (^5]) . 

By Skohorod representation Theorem |Bil99j . and since C{[Q,T]; H~^{T'^)) is 
separable, there exists a new probability space (f2,-F, P) and some random va- 
riables 



23 



with respective law and P such that p in C([0, T]; i7-''(T'*)), P a.s. 

Let (f S C'^(L^(T'*)) be regularizing and subquadratic according to ([M)) . By 
Corollary [9] and the L^-bound of Proposition [TOl we have for < < • • • < 



E 



V{p'{t))~^{p^{s)) 



< Ce 
(66) 

with a constant C depending on the constant C* in on Cip, on sup^^g ||/o Hl^ 
and on the supremum of ij}. Since p'^ and p"^ have the same law, this is still true 
if p"^ is replaced by . Assume furthermore that Lp is bounded and continuous 
from iJ^''(T'^) into M, then it is easy to take the limit £ — >■ in ([55)) and to 
obtain 



V{p{t)) - - / ^^(p(a))da 



V'(p(si),...,p(s„)) =0. (67) 



E 



The additional hypothesis on ip can be relaxed. Indeed, thanks to Proposition 
1101 we can approximate every subquadratic and regularizing functions by func- 
tions in ChiH'"^ {T'^)) which are subquadratic and regularizing with a uniform 
constant in p4p and which converge pointwise. 

We have thus proved that P solves the martingale problem associated to ^ 
with subquadratic and regularizing test functions. In particular, for all such ip: 

M^{t) = ip{p{t)) - [ ^^{p{s))ds, t > 0, (68) 
Jo 

is a martingale with respect to the filtration J-g generated by (p(s)). The 
quadratic variation of is (c/ (|33p ) 

Let D(p{p) (E) Dfip) denote the bilinear form 

{h,k)^{h,D^{p))ik,Dipip)) 
on L^{T'^). By (113), we have 



ifv''(p)-2^(p(s))Jfv'(p(s)) 

We deduce that 

M{t) = p{t) - p(0) 



^Tmce[pQ^/^)D^{p) ^ D(p{p)pQ^/^) 



pQ^'^D^ip) 



L2 



-^Pis) + -Fp{s)ds, t > 0, 



is a martingale with quadratic variation j p{s)Q^^^ yPis)Q^^^j ds. Thanks 

to martingale representation results (see for instance |DPZ92| ). up to a change 
of probability space, there exist a cylindrical Wiener process W such that 



Pit) 



p(0) - / Apis) + lFp{s)ds ^ f p{s)Q^^^dW{s), t > 0. 
Jo 2 Jo 
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It is well known that this equation has a unique solution with paths in the 
space C{[0,T]; H~'^{M.'^)). This can be shown for instance by energy estimates 
using Ito formula after a suitable regularization argument. Moreover pathwise 
uniqueness implies uniqueness in law and we deduce that P is the law of this 
solution and is uniquely determined. Finally, by uniqueness of the limit, the 
whole sequence {P'^) converges to P weakly in the space of probability measures 
on C([0,r];if-"(R'^)). 
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